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Abstract 

The Wigner function for one and two-mode quantum systems is explicitely expressed in 
terms of the marginal distribution for the generic linearly transformed quadratures. Then, 
also the density operator of those systems is written in terms of the marginal distribution 
of these quadratures. Some examples to apply this formalism, and a reduction to the usual 
optical homodyne tomography are considered. 
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1 Introduction 



The homodyne measurement of an electromagnetic field gives all the possible combinations of the 
field quadratures belonging to the group 0(2), by just varying the phase of the local oscillator. 
The average of the random outcomes of the measurement, at a given local oscillator phase, is con- 
nected with the marginal distribution of the Wigner function, or any other quasi-probability used 
in quantum optics. In the work |I[ it was shown, indeed, that the rotated quadrature distribution 
may be expressed in terms of the Wigner function || as well as in terms of the Husimi Q-function 
H and the Glauber-Sudarshan P-representation [[|], p. The result of Ref. [I] was based on the 
observation that relations of the density matrix in different representations to the characteristic 
function, intended as the mean value of the displacement operator creating a coherent state from 
the vacuum [|J, ||, may be rewritten as relations where the characteristic function is the mean 
value of the rotated quadrature phase. It gave the possibility to express the Wigner function in 
terms of the marginal distribution of homodyne outcomes through the tomographic formula. The 
essential point of the obtained formula is that the homodyne output distribution may be measured 
and the corresponding inverse Radon transform produces the Wigner distribution function, which 
is refered as the optical homodyne tomography J7|. The density matrix elements, in some repre- 
sentations, can also be obtained by avoiding the Wigner function and then the Radon transform 

I, 1, 0- 



The tomographic formula of Ref. [1] has been presented in an invariant form in Ref. |K|. In 
that work the system density operator was expressed as convolution of the marginal distribution 
of the homodyne output and a kernel operator. In Ref. [IT| the tomographic formula of Ref. [1] 
was extended to express the Wigner function in terms of the marginal distribution of the gener- 
alized quadrature obtained through canonical transforms belonging to the Lie group ISp(2,R). 
Essentially this generalized quadrature differs from the one used in the standard tomography since 
it depends upon three real parameters instead of only one (really, a scaling transformation is also 
present in the homodyne measurement because position and momentum have different units, but 
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the scaling cannot be independently controlled of the rotation). 

The aim of this paper is to obtain the invariant form of the density operator of the quantum 



system in terms of the discussed [11] marginal distributions for the generalized quadratures, which 



is analogous to Ref. fL0| , and to extend the analysis to the multimode case. 



Recently the quantum tomography was also considered to investigate material systems |L2 
besides the optical ones, thus our alternative approach could be useful for a better characterization 
of their quantum state since it allows to get overcomplete information. 

The structure of the article is the following. In Section 2 we derive the density operator as a 
convolution of the marginal distribution for the generic linear quadrature and a kernel operator, 
and discuss the properties of the obtained formula in number state basis, coherent states basis 
and coordinate representation. In Section 3 we consider the obtained formula applied to some 
examples of interest, such as the ground state of a harmonic oscillator, the Schrodinger-cat state, 
and the thermal equilibrium state of an oscillator. In Section 4 we contrast the obtained formula 
with the tomographic one for the density operator of Ref. [l(| . In Sections 5 and 6 the approach 



will be extended to the multimode case. Some applications of the developped formalism as well 
as practical implementation and further developments will be discussed in Sections 7 and 8. 

2 Density Matrix and Marginal Distribution 

Let us consider the quadrature observable X, which is a generic linear form in position q and 
momentum p 

X = fxq + up + S (1) 

with n,v,8 real parameters (their physical meaning will be discuss later on), then it is possible 
to get the density matrix elements from the marginal distribution avoiding the evaluation of 
the Wigner function as an intermediate step. For this pourpose we start from a well known || 
representation of the density operator 

P = I —W(a)f(a) (2) 

J 71 



where the Wigner function W(ct) is a weight function for the expansion of the density operator in 
terms of the operator T(a) which is defined as the complex Fourier transform of the displacement 
operator D 

f(a)= /^L>(£) eX pK*-a*0- (3) 

J 7T 



Following the lines of Ref. |TT[ we may write the marginal distribution w for the generic quadrature 
of Eq. (|) as 

v>{X,n,v,S) = J e-^ x -^-*) W ( q , P )^^, (4) 

where the canonical coordinate q and p are related to a by 

q + ip 

Eq. (^) shows that w is function of the difference X — 5 = x, so that it can be rewritten as 

w(x,n,u) = J e - ik ^-^W(q,p)^^; (6) 

and by means of the Fourier transform of the function w one can obtain the relation 

W(q,p) = (2it) 2 z 2 w(z, -zq, -zp) (7) 

where — zq, —zp, z are the conjugate variable to \l, v, x respectively. It is worth to remark that 
in this case the connection between the Wigner function and the marginal distribution is simply 
guaranteed by means of Fourier transform instead of the Radon one. Then inserting Eq. (|7|) into 
Eq. (0), expressing the marginal distribution in terms of the Fourier transform, we have: 

p = J ^^f(q,p) J dxdfidv z 2 w(x,/j,,iy)e- izx+i ^ zq+iuzp , (8) 

or, in a compact form 

p = J dxdjidu w(x, p, v)K^ v , (9) 
where the kernel operator K^^ is given by: 

k w = j ^ 2 z 2 ^ Zq+lUZp - lZX T{q,p) 

1 . z 2 e - izx e -^ y - i ^ e ^ u+i ^ e -4^ 2+v2 ) (10) 
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where we have used Eq. (|3]) and the Baker-Hausdorff formula. The fact that K^ v depends 
on the z variable as well (i.e. each Fourier component gives a selfconsistent kernel) shows the 
overcompleteness of information achievable by measuring the observable of Eq. ([!]). Thus, Eq. 
(H) can be useful to completely determine the properties of the system, i. e. the density operator, 
from the probability distribution of the experimental data; outcomes of a set of measurements 



such as described in Ref. 11 



Let us now consider the expression of the kernel operator (|T0|) in various representations. First, 
in the coherent states basis we have 

(a|iT 1/3) = iA-e^ (v - v)a ^ ( ^e-T(^^e-^4+«^ (11) 
while in the number states representation we get 

(n + d\K^\n) = * V^e-^ 2+ ^ 

X (n-l)\ l\ {l + d)\ ' { } 

and finally in the coordinate representation we obtain 

W\Kuu\q") = —z 2 e- izx e lz2 ^ /2 e izm '5{q" - zv - q'). (13) 
2n 

If, and only if, the kernel operator is bounded every moment of the kernel is bounded for all 
possible distributions w(x, /i, v\ then, according to the central limit theorem, the matrix elements 
of Eq. (^) can be sampled on a sufficiently large sets of data only in the number or coherent states 
basis, as can be evicted from the expressions (0), (|T2l) and (|i~3|) . That is in perfect agreement 



with the results of Ref. 10 



3 Examples 

In this Section we will consider some examples of interest to check the validity of Eq. (Q). We 
first consider the case of the harmonic oscillator's ground state for which we have [Q : 

w(x, fi, v) = [ir(ji + v)Y xl2 exp[ — T — —}. (14) 
5 



In the coherent states basis of Eq. (|TT| ) we obtain: 

\a\ 2 \f}\ 2 

dxd\idv w(x, fi,i , )(a\K^ u \P) = e 2 2 , (15) 

which exactly is the product (a\0)(0\(3) as we expected. As another example we consider the 
Schrddinger cat state |13j of the type discussed in Ref. [0 

| a + ib) + | a — ib) 



{2[1 + cos(2a6)exp(-26 2 )]} 1 /2 



(16) 



with a and b arbitrary real numbers, for which we have a marginal distribution given by |TT 

(x — fia) 2 + b 2 v 2 



w(x, fl, v) 



(T 


r 1 " 


1/2 1 


V 7T y 


-/I 2 + V 2 - 


[1 + cos(2a6) exp(-26 2 )] 



cxp 



X 



jcosh 



4z/6(x — fia) 



fi 2 + v 2 



+ cos 



2&(2/ia; - a(/i 2 - v 2 )) 
fi 2 + v 2 



]}■ 



fi 2 + V 2 



(17) 



Then, from Eq. @ we have in the coherent states basis of Eq. (]TTf) 



J dxdfidv w(x, n,v){a.\K^, 



-(a 2 +b 2 )-(|a| 2 +|/3| 2 )/2 



3 (a*+/3)(a+ib) _j_ e a* (a+ib)+f3(a-ib) _|_ g a* (a- ib)+/3(a+*) _|_ e (a*+/3)(a-i&) j /-^ 



[1 + cos(2a6)e- 262 ] 

which exactly is the product as we expect. Finally, we consider the application of Eq. 

(|^) to the thermal equilibrium state of the oscillator at temperature T. In this case the Wigner 
function is given by [14{] : 



W T {q,p) = 2 tanh (|^) exp [ - {q 2 + p 2 ) tanh (^)]- 



(19) 



From this equation the marginal distribution can be easily obtained as follows [|Tl 

w(x, fi, v) 



1 1 

2"7T \i 



/X V 
dp W T { p,p) 
u u 



A 



,1/2 



7r(/X 2 + Z/ 2 )" 



/i /X 

exp [—A 



x 



fi 2 + z/ 2 



(20) 



where we have introduced A = tanh {^). Now, we use Eqs. ( PD and ( |20|) to calculate 



>2fcT^ 



/ dxdfidv w(x, fi,u)(a\K i 



fl,V I 



A 



1 - A 



2 h — r] ex Ph — r"*/? - 
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a 



101 s 



(21) 



This result is the same that we could have by directly performing (a\px\P) with the thermal 
density operator p? given by jl4| 



p T = (1 - e -^/ feT ) J2 \n) (n\e- nhw/kT . (22) 

n 

4 Comparision with the usual tomographic technique 

A relation betweeen the density operator and the marginal distribution analogous to that of Eq. 
(|]) can be derived starting from another operator identity such as pi 

d 2 a r 



P 



f —Tr{pD(a)}D~\a) (23) 

J 7T 

which, by the change of variables p = —y2^sma, v = \^2^R,ea, becomes 

P = — f dpdv Ti{pe- l± }e l± = — f dpdv Tr{pe- i£ }e i& (24) 
2% J 2n J 

where x = X — 5. The trace can be now evaluated using the complete set of eigenvectors {\x)} of 

the operator x, obtaining 

Tr{p e - i& } = Jdx w(x,p,u)e- ix (25) 

then, putting this one into Eq. (p4|), we have a relation of the same form of Eq. @ with the 
kernel given by 

£ = — e ~ ix e ix = —e-^e^+^P^ (26) 
2tt 2tt 

which is the same of Eq. (|I0|) setting z = 1. It means that we now have only one particular Fourier 
component due to the particular change of variables (the most general should be zp = —\/2^sma 
and zv = \/2^R.ea). 
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In order to reconstruct the usual tomographic formula for the homodyne detection |1Q 
need to pass in polar variables, i.e. p = — rcos0, v = — rsin0, then 

x — > — rx^ = — r[q cos (j) + psin0]. (27) 

Furthermore, denoting with x$ the eigenvalues of the quadrature operator x^, we have 

Ti{p e - ix } = Tr{pe irx *} = J d X4> w(x^<P)e irx * (28) 



and thus, from Eq. ( |24| ) 

p = J dcpdx^ w{x^4>)K<t, (29) 

with 

Ka ) = — [ dr re ir ^- £ ^ (30) 
Itx J 

which is the same of Ref . |10| . Substantially, the kernel of Eq. (|30| ) is given by the radial integral 
of the kernel of Eq. (|26|), and this is due to the fact that we go from a general transformation, 
with two free parameters, to a particular transformation (homodyne rotation) with only one free 
parameter and then we need to integrate over the other one. 

5 Two- Mode Tomography 

In this Section we discuss an extension of the above approach to multimode systems and, for 
simplicity and applications, we concentrate on the two-mode case. From now on we use the 
convention that the vector symbol represents only two dimensional vectors. 

5.1 Two-Mode Quasiprobability and Marginal Distribution for Two 
Squeezed and Shifted Quadratures 

We introduce the vector operators q = (qi, q 2 ), V = {P11P2) and then the following observables 

= A (q,$) T + (5,5') (31) 
in which A is a real symplectic 4x4 matrix 



AoA 



a: a 



( 10^ 
1 
-10 

V -100/ 



(32) 



and (5,5') = (Si, 62, 6[, 8' 2 ) is a real c-number four-vector corresponding to the shifts of the 
four quadratures. Furthermore, the components Xi, X2, Yl, Y2 are related to the inhomogeneous 
symplectic group ISp(A,R). Then we introduce the marginal distribution for the observable 
X = (Xi,X 2 ) due to relations which, obviously, are two-mode generalizations of the connection 



between the characteristic function and the marginal distribution, given in Ref. || for one-mode 
case. We have firstly the characteristic function as 

x (k) = ( e ^) = Tr{pe 4 ^}; kk = hX, + k 2 X 2 (33) 

where k — (ki, k 2 ), and p is the density operator of the two-mode system. In terms of the Wigner 
function W(q,p) of the system it is 

X(k) = J e^W(q } p) dqdp. (34) 

Here X is, in principle, any observable, but to be concrete below we will consider this vector to 
have components 

Xi = pq + vp + Si 

X 2 = pq + vp + S 2 (35) 

with p = (Au,Ai 2 ); V = (Ai 3 ,A 14 ); p' = (A 21 ,A 22 ); V = (A 23 ,A 24 ). The marginal distribu- 
tion w(X, p, i>, p', P 1 , S) of the vector variable X, depending on the parameters of the symplectic 
transformation A and shift vector 5, is given by the relation 

w(X,p,v,p',i?,5) = ^jdk X (k)e- iM 

= J dkdqdpWi&pjexpl-ik^Xx -pq-vp-b x ) 

- ik 2 (X 2 - pq- up- 5 2 )]. (36) 

From this formula it is clear that also in the two-mode case the marginal distribution depends 
on the difference X — 5 = x, so that the replacement w(X , p, v, p! ', V ', 6) — > w(x, p, V, p', V) is 
possible, with the normalization 

w(x, p, P, p', if) dx = 1. (37) 

Our aim is to express the Wigner function W(q, p) in terms of the marginal distribution probability 
w(x, p, v, p! , z/'), i.e. to invert the formula (^). For this purpouse we perform a Fourier transforms 
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of Eq. (|36|) and, after some minor algebra [TTJ , we arrive at 



Wi-m/z!, -n/z x ) = (2Txfz i 1 e- i{ f 1 ' MH)z2,Zl w{z, m, ft, V) (38) 

where the vectors m = (mi, m.2); ft = (ni, n 2 ); z — (zi, z 2 ) are the conjugate variables to jl, V and 
x respectively. 

5.2 Two-Mode Quasiprobability and Marginal Distribution for One 
Squeezed and Shifted Quadrature 

We will now discuss the connection between the Wigner function W(q,p) and the marginal dis- 
tribution w(xi,fl, V) for only one quadrature X 1 . This latter is related to the previous by 

w(xi,jl, v) = Jw(xi 1 X2,[I,i',fl',if)dx2- (39) 

Inserting in this one the expression of w given in Eq. (|36D, we may see that the marginal distri- 
bution w is connected with the characteristic function = x(.ki,k 2 = 0); where x(.k\,k 2 ) is 
given in Eq. ([!3|). Repeating step by step the procedure that leads to Eq. ( pB]) we have in this 
case 

W(-m/zi, -n/zx) = (27r) 3 z^w(z 1 ,m,n). (40) 



This formula is similar to the one-mode case discussed in Ref. |TT[ and in the previous section, Eq. 
(0). The possibility to relate the Wigner function of the two-mode system with different marginal 
distributions, w(x) or w(xi), is connected with the use of different quantities of information ob- 
tainable from measurements. Of course, the two-dimensional distribution function contains much 
more experimental information than the one-dimensional. However, in both cases, information are 
overcomplete to determine the Wigner function of the quantum state, but must be self consistent. 
In concrete situations these varieties of descriptions might give a freedom to choose the convenient 
type of measurements. 



10 



6 Invariant Expression for Two-Mode Density Operator 

We start from a generalization of Eq. (0) to the two-mode case 

d 2 a\d 2 a2 



n 2 



W(a)T(a) (41) 



with a = (01,0^2) and ctj = (qj + ipj) /V%] j = 1,2. Then inserting in this one the Wigner function 
of Eq. (|38|), and expressing the marginal distribution in terms of its Fourier transform, we have 

p = J ^^^(q,p) J dxdpdPzfe iZ2 ^' ?+ ^w(x,p,P,p,i/)e iz ^ +iZli; P- igs (42) 

that can be written as 

p = J dxdpdvw(x, JX, V, p, v)K^^i^ (43) 
with the kernel operator explicitely given by 

f> 1 4 — izx 

ix,v,pl,V — ( 27r )2 1 

x exp|-^=[zi(z/ + ip) + z 2 {v + ijH)]a— -^=[zi(i> — ip) + z 2 {v - ijH)]a} (44) 

with a = (01,02) the vector operator describing the two modes. In the case where the Wigner 
function is given by the marginal distribution of only one quadrature, as in Eq. (fRJ), the same 
procedure leads to an expression similar to Eq. (|43|), but with a different kernel operator 

Kp,p= j2^4^ izixl eiq>{^[(u + ip)S- (v-ip)a]} (45) 

which is very similar to that of Eq. (|T0|) . It should be noticed that as a direct extension of the 
arguments of Sec. 2, also the kernels in Eqs. fl44|) and ( fl5|) are not bounded in the quadrature 
(or coordinate) representation, then in this case it is not possible to sample the density matrix 
elements. 

7 Applications 



Here we shall consider some applications of the presented tomographic scheme. We are aware 
that the crucial point might be the practical achievement of the generic linear quadrature such 
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as in Eq. (pq). To consider, however, a measurement scheme as an optical implementation of the 
developped formalism, let us go back, at first, to the one dimensional case. The quadrature of 
Eq. ([]]) could be experimentally accessible by using for example the squeezing pre-amplification 
(pre-attenuation) of a field mode which is going to be measured (a similar method in different 



context was discussed in Ref. |15|). In fact, let a be the signal field mode to be detected, when it 



passes through a squeezer it becomes a s = a cosh s — a^e ld sinh s, where s and 9 characterize the 



complex squeezing parameter ( = se jig] . Then, if we subsequently detect the field by using the 
balanced homodyne scheme, we get an output signal proportiopnal to the average of the following 
quadrature 

%,)= * (a s e-^ + at e ^), (46) 

where <fi is the local oscillator phase. When this phase is locked to that of the squeezer, such that 
<f) = 6/2, Eq. (H) becomes 

E(</>) = (ae^ /2 [cosh s - sinh s] + aV^cosh s - sinh s]) , (47) 

which, essentially, coincides with Eq. ([I]), if one recognizes the independent parameters 

/i = [coshs — sinhs] cos(6 l /2); v = [coshs — sinhs] sin(#/2) . (48) 

The shift parameter S has not a real physical meaning, since it causes only a displacement of the 
distribution along the X line without changing its shape, as can be evicted from Eqs. (|) and (|6|). 
So, in a practical situation it can be omitted. 

In the two-mode case it could be interesting to use the connection between the Wigner func- 
tion and the marginal distribution of only one quadrature in the case of heterodyne detection 
(particulary used to detect multimode squeezed states). We may refer to this scheme as the op- 
tical heterodyne tomography. In fact, in balanced heterodyne detection the measured current is 
determined by [16| 



E(4>) = 4=[e**oI + e-^ai] + %e**o5 + e~^a 2 ] . (49) 
v 2 v 2 
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If one uses squeezed pre-amplification (pre-attenuation) for each single mode, as discussed 
above, both E% and E2 are indipendent variables, however, we now have only one phase, the 
local oscillator phase 0; thus, in order to have a fourth independent variable which would allow a 
complete reconstruction of two-mode quasi-probability, we could consider the phase shifter used 
in the first step of the measurement procedure. Since the two modes have different frequencies, 
the phase shifter induces different phase changes in the two modes; for example 



<2i — > a x e %dx ; d 2 — > a 2 e l6 ' 2 ; 



(50) 



where the difference Q\ — 62 depends on the length of the optical path in the phase shifter. Upon 
using this requirement, the measured quadrature ( [Hp will become 



V2 



E2 
V2 l 



(51) 



that corresponds to the quadrature X% of Eq. ( |3"5l ) with \x = (E\ cos(0 + #1), E 2 cos(0 + 62)) and 
V = (Ex sin(0 + 81), E 2 sin(0 + ^2))- For the shift parameter Si the previous considerations, made 
in one mode case, still hold. 

We would now present some examples of the probability w measurable with the above scheme. 
We first consider the most general two-mode squeezed states described by the Wigner funtion of 



the form 17 



W(q,p) = (detM)- 1 /2exp[-i((g,p) - ((£ ^M' 1 ((g, p) - Umf] (52) 
where the real symmetric dispersion matrix A4 has ten variances 

^c c i3 = ^((q,P) a {q,p)p + {q,P)p{q,P) a ) - ({q,P)a)({q,P)p) «, P = 1,2,3,4. (53) 

By integrating only on k\ with ft2 = 0, we obtain the marginal distribution of only one quadrature 
Xi as given in Eq. ( p6|) 

1 



w{xi,n, v) 



dk\ I dqdpW(q,p)e 



(2tt) 3 

[2^iT)M{fl^) T r l/2 



cxp 



-iki(xi—flq—up) 



.1 1 



2(fl,u)M(fi,u) 



(54) 
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where we have assumed ((q,p) a ) = for all a. This distribution is a Gaussian for the quadrature 
variable. 

As a second example, we consider the two-mode Schrodinger-cat state |1| defined as 



\A±) = N±(\A)±\-A)) 



(55) 



with 



,\A\*/2 



a |A|»/2 



2\/cosh \A\ 2 



2\/sinh \A\ 2 



(56) 



The Wigner function of this state is 



W A± (q,p)=Nl [ W {A ^ =A) (q } p)±W [A§= _ A) (lp) 



± W { _ XiS=X) (q,p)±W l 



(-A,B=-A) 



(q,p)] 



where 



W ab = 4ex P 



-2aa* + 2Aa* + 2B*a - AB* 



\A\ 2 \B\ 



(57) 



(5f 



is the Wigner function for two-mode coherent state |H| with 

q + ip 



a 



V2 



(59) 



Let us now consider the even coherent state, then to get the marginal distribution, we use, as in 
the previous case, Eq. ([36]) obtaining 



x jexp 



2n- 1 / 2 N 2 + { -x 2 - {v x P x + /ixQO 2 - {v 2 P 2 + /i 2 Q 2 ^ 2 
exp 



jl\ 2 + \u\ 2 



\fl\ 2 + |z7| 2 

(P 2 + Ql){v 2 + - (P| + Q 2 )(u 2 + tf) + 2(// 1 P 1 - u 1 Q 1 )( f x 2 P 2 - u 2 Q 2 ) 



\v>\ + M 



x cos 



+ exp 



2(/iiPi + /i 2 P 2 - f 1Q1 - ^2^2)^1 



|/z| 2 + |z?| 2 
-2(z/ 1 P 1 + /iiQi)(z/ 2 P 2 + // 2 Q 2 ) 



I A* I + M 



cosh 



2(i/i Pi + V2P2 + H1Q1 + ViQ-i)^ 



w\ + \ v \ 



}• (60) 



where we have set 



.4 



Q + iP 
V2 ' 



(61) 
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In the case of \ii = v<i = and Q2 = P2 = 0, the marginal distribution of Eq. (|60| ) reduces to the 
one-dimensional marginal distribution of Eq. (32) of Ref. (Tl| 0. 

8 Conclusions 

We have shown that the optical tomography may be extended to two-mode case, and the new 
observables introduced, which are generic linear forms in quadratures, give the possibility of using 
various measurements, different from the usual ones, to determine the Wigner function of the 
system under study in terms of the marginal distributions. To better understand the meaning 
of the measure of the observable X of Eq. (Q), we recall that this symplectic transformation 
could be represented as a composition of shift, rotation and squeezing. Along this line we have 
presented some schemes which could be implemented in optics. To be more precise, the shift 
parameter does not play a real physical role in the measurement process, it has been introduced 
for formal completeness and it expresses the possibility to achieve the desired marginal distribution 
by performing the measurements in an ensemble of frames which are each other shifted; (related 
method was early discussed in Ref. p0| ). In an electro-optical system this only means to have the 
freedom of using different photocurrent scales in which the zero is shifted by a known amount. 

The tomographic formalism is presented here in an invariant form for the density operator 
and it may be suitable for further numerical analysis. It is worth to remark that all the formulae 
of two-mode case may be obviously rewritten for N-mode case with arbitrary N. The presented 
extension gives the possibility to include the procedure of heterodyne detection (two-mode case) 
as well. Furthermore, the extension of the tomographic technique to multimode systems could 
also be useful to obtain the (complete) quantum information about a system in an indirect way 
21 . 

Another conclusive remark is related to the group structure of the presented tomographic ex- 
tension. Using generic observables which are linear in quadratures, we apply a symplectic group 

transform to an initial quadrature component. In one mode optical tomography this transform 
2 Up to a misprint in this reference. 
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belongs to the rotation subgroup of the symplectic group ISp(2,R). Thus the presented con- 
struction of the Wigner function (density operator) may be reformulated as the group problem 
for a particular symplectic group and we could call our extension as " symplectic optical tomogra- 
phy" . On the other hand the scheme might be generalized to other Lie groups different from the 
symplectic one. We will develop these points in future papers. 
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